In Search of Circumstellar Disks Around Young Massive Stars 



Luis A. Zapata^'^, Luis F. Rodriguez^, Paul T. P. Ho^'^, Henrik Beuther^ and Qizhou 

Zhang^ 

IzapataOastrosmo . unam . mx 

Received ; accepted 



Submitted to The Astronomical Journal 



^Centre de Radioastronomi'a y Astrofisica, UNAM, Apdo. Postal 3-72 (Xangari), 58089 
Morelia, Michoacan, Mexico 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, 
USA 

■^Academia Sinica Institute of Astronomy and Astrophysics, Taipei, Taiwan. 



- 2 - 



ABSTRACT 

We present 7 mm, 1.3 cm and 3.6 cm continuum observations made with 
the Very Large Array toward a sample of ten luminous IRAS sources that are 
believed to be regions of massive star formation. We detect compact 7 mm 
emission in four of these objects: IRAS 18089-1732(1), IRAS 18182-1433, IRAS 
18264-1152 and IRAS 18308-0841 and for the first time find that these IRAS 
sources are associated with double or triple radio sources separated by a few 
arcseconds. We discuss the characteristics of these sources based mostly on their 
spectral indices and find that their nature is diverse. Some features indicate that 
the 7 mm emission is dominated by dust from disks or envelopes. Toward other 
components the 7 mm emission appears to be dominated by free-free radiation, 
both from ionized outflows or from optically thick H II regions. Furthermore, 
there is evidence of synchrotron contamination in some of these sources. Finally, 
we found that the sources associated with ionized outflows, or thermal jets are 
correlated with CH3OH masers. The precise determination of the nature of 
these objects requires additional multifrequency observations at high angular 
resolution. The 3.6 cm continuum observations also revealed seven UCHII regions 
in the vicinity of the sources IRAS 18089-1732(1) and two more in the source 
IRAS 18182-1433. We show that the small photoionized nebulae of these UCHII 
regions are produced by early B-type stars. 

Subject headings: stars: massive star formation - stars: pre-main sequence - ISM: 
radio thermal jets - Individual: IRAS 18089-1732 - Individual: IRAS 18182-1433 - 
Individual: IRAS 18264-1152 - Individual: IRAS 18308-0841 
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1. Introduction 

Massive stars (M>10 M0) are fundamental in determining the physical and chemical 
evolution of the galaxies. They produce a dominant fraction of the heavy elements, generate 
large amounts of UV radiation during their lives, and inject turbulent energy into the 
galactic ISM. However, the formation of massive stars is a poorly understood phenomenon. 
Schematically, there are two scenarios that have been proposed to explain the formation of 
massive stars. The first is the accretion scenario, (Garay & Lizano 1998, Norberg & Maeder 
2000, and McKee & Tan 2002) similar to that operating for low mass, where a dynamical 
collapse of molecular clumps, results in the formation of disks and molecular outflows, 
which leads to the formation of a stellar core. In contrast to this, the coalescence scenario 
(Bonnell, Bate & Zinnecker 1998; Stabler, Palla & Ho 2000; and Bally & Zinnecker 2005) 
proposes that massive stars form by the merging of stars of low and intermediate masses. 
The possible existence of this coalescence scenario receives support from two facts: i) ior a, 
spherically symmetric geometry, the radiation pressure can reverse the infall process once 
the central star reaches about 10 Mq (Wolfire & CassineUi 1987; see, however, Yorke & 
Sonnhalter 2002 for models that can produce stars with larger masses) and ii) most massive 
stars are born in the center of dense clusters of low and intermediate mass stars (Bonnell 
et al. 1998, and references therein). While the accretion scenario implies the existence of 
well formed disks and jets, this is not the case for in the coalescence scenario, where disks 
and jets are expected to be disrupted during the merging. 

In recent years there has been a large effort to find observational evidence to 
discriminate between the two scenarios. Searches have been carried out to find evidence 
for rotating disks using VLA ammonia observations. In two out of twelve objects studied, 
both B-type protostars, Zhang, Hunter & Sridharan (1998) and Zhang et al. (2002) found 
evidence for signatures of rotating disks. Pestalozzi et al. (2004) have interpreted VLBI 
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observations of methanol masers in NGC7538 IRSl in terms of an edge-on Keplerian 
disk extending to a radius of ~1000 AU and orbiting a 30 Mq protostar. However, the 
methanol masers seem to be often associated with shocked molecular gas (De Buizer 2003; 
Kurtz 2004). Gibb et al. (2004) have tried imaging dusty disks using subarcsecond BIMA 
millimeter observations (at 2.7 mm and 1.4 mm wavelengths) toward eight massive objects. 
In most of the objects they found that even at 2.7 mm the dominant emission mechanism 
is free-free continuum from a thermal jet. Moreover, if other OB-type stars have been 
formed previously in the region, the object studied may be in the vicinity of bright free-free 
emission that will make it very difficult to search for the relatively weak mm emission from 
the disk. Finally, the shorter photoevaporation timescale of disks, the large distances (> 
Kpc), and the heavy extinctions (Ay > 100) have also been inherently difficult problems in 
the search for disks around high-mass protostars. 

At present, there is a hst of early B-type protostars that have been associated with 
possibles circumstellar disks, some of them are: the BN object (Jiang et al. 2005), Cepheus 
A HW2 (Patel et al. 2005), IRAS 20126+4140 (Cesaroni et al. 1997; 1999; Zhang, Hunter & 
Sridharan 1998; Zhang et al. 2000), IRAS 18089-1732 (Beuther et al. 2004b, 2005b), IRAS 
18182-1433 (Beuther et al. 2006; Zhang et al. 2006, both in preparation.), G192.16-3.82 
(Shepherd, Claussen, & Kurtz 2001); W33A, AFGL 2591 & NGC 7538 IRS9 (Van der Tak 
& Menten 2005), AFGL 5142 (Zhang et al. 2002), AFGL 490 (Schreyer et al. 2002 and 
2005 in preparation). 

We note that sources IRAS 18089-1732 and IRAS 18182-1433 are remarkable examples 
of massive protostars that are associated with possible circumstellar disks that will be 
further discussed here. 

In this paper, we present a centimeter and millimeter wavelength continuum study of 
ten IRAS sources that are thought to be regions of massive star formation. The purpose of 
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this research is to identify the nature of these massive sources and to select a subset of good 
candidates for having associated a circumstellar disk. 

2. The Sample 

The ten IRAS regions observed were selected from the sample of 69 high mass 
protostellar objects (HMPOs) of Sridharan et al. (2002) using the following criteria: high 
luminosity (~2xl0^ Lq) suggesting the presence of OB-type protostars; compact size in 
their mm emission (< 20"x20"); proximity to the Sun (a few Kpc or less); bright flux 
density at 1.2 mm (>200 mJy), and weak centimeter emission (<1 mJy). The criterion of 
weak centimeter emission was included to search for sources where the continuum emission 
at 7 mm and shorter wavelengths from possible disks is less likely to be contaminated by 
free-free emission. In Table 1 we show the main properties of the ten objects selected. 

3. Observations 

The observations were made with the Very Large Array of NRAO^ in the continuum 
mode at 7 mm during 2004 April 1, and at 1.3 cm and 3.6 cm on 2004 May 17. At these 
epochs, the VLA was in its C configuration. The absolute flux calibrator was 1331+305. In 

Tables 2, 3 and 4 we present a summary of the observations that will be discussed here, 
including the measured flux density of the phase calibrators. 

As mentioned above, during the first epoch, wc made 7 mm continuum observations 
toward the ten selected massive young regions. We detected continuum emission at or above 

^The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 
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a 3-(7 level of ~0.5 mJy in four of the regions: IRAS 18089-1732(1), IRAS 18182-1433, 
IRAS 18264-1152, and IRAS 18308-0841. During the second epoch, we observed at 1.3 and 
3.6 cm wavelengths toward these four regions. 

The data were edited and calibrated in the standard manner using the software package 
Astronomical Image Processing System (AIPS) of NRAO. Clean maps were obtained using 
the task IMAGR of AIPS. For the 7 mm continuum maps we used the ROBUST parameter 
set to 5, that corresponds to natural weighting, to achieve maximum sensitivity in each 
continuum image. In the 1.3 and 3.6 cm continuum maps we used the ROBUST parameter 
set to 0, for optimal compromise between sensitivity and angular resolution. The contour 
and grey scale maps at 0.7, 1.3 and 3.6 cm of the four regions are shown in Figures 1, 2, 4, 
5, 7, 8, 10, and 11. The resulting rms noises and synthesized beam parameters for the 0.7, 
1.3 and 3.6 cm continuum images are also given in Tables 2, 3 and 4. 

Additionally, in all 7 mm maps we have also applied a uv-tapering between 100 to 200 
kA in order to increase the signal-to-noise ratio of the extended emission (> 1") in these 
regions. 

4. Results and Discussion. 

In what follows, we discuss each of the four sources detected at 7 mm separately. Our 
discussion of the spectral indices is based on the assumption that the flux densities did not 
change between the 2004 April 1 and 2004 May 17 observations. 
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4.1. IRAS 18089-1732(1) 

This source is previously known as IRAS 18089-1732, however we have added a 
numbering (1) at the end of the name to distinguish between this component and 
component IRAS 18089-1732(4) 

This region contains a strong and extended (~ 11") 1.2 mm continuum source, H2O 
and Class II CH3OH maser spots, and a compact (~ 1") 3.6 cm continuum source that is 
centered at the middle of the extended 1.2 mm continuum source (Beuther et al. 2002a, 
2002c). In recent millimeter and submillimeter observations made with the Submillimter 
Array (SMA), Beuther et al. (2004b, 2005b) found a 1.2 mm continuum compact core 
elongated in the southeast-northwest direction, while in the 860 /im continuum observations 
they detected emission more compact, that is coincident with the 3.6 cm source and the H2O 
maser (see Beuther 2004b, Figure 1). Furthermore, Beuther et al. (2004b, 2005b) found a 
molecular line forest associated with this source, and particularly the HCOOCH3(20-19) line 
appears to show evidence for the signatures of rotation. Finally, IRAM 30 m single-dish SiO 
J=2— )-l (Beuther, priv. comm.), and SMA SiO J=5^4 (Beuther et al. 2004b) observations 
show evidence of outflowing molecular gas emanating from this region. 

The 7 mm maps (see Figures 1 and 2) show that this source is actually a double 
source with components separated approximately by ~ 2" (7200 AU, at a distance of 3.6 
Kpc) in the northeast and southwest direction. We refer to these components as IRAS 
18089-1732(l)a and IRAS 18089-1732(l)b. 

In Figure 2 we have overlayed the 1.2 and 0.86 mm continuum maps from Beuther 
et al. (2004b) with our 7 mm continuum maps. Figures 2a and 2b show that the 1.2 and 
0.87 mm continuum cores found by Beuther et al. are associated with the compact 7 mm 
continuum source IRAS 18089-1732(l)a. We now discuss each subcomponent separately. 
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4.1.1. IRAS 18089-n32(l)a 

The radio source IRAS 18089- 1732(1) a is resolved in the 7 mm maps. It has 
deconvolved dimensions of 2'!?, ± OM x V!\ ± 0^'4 and a P.A. = 28° ± 15°. In the 1.3 cm 
continuum maps we detect the counterpart of IRAS 18089-1732(l)a (see Figure lb). This 
source is also resolved and has deconvolved dimensions of 1''4 ± 0''2 x 0''9 ± 0''2 and a 
RA. = 32° ± 13°. The position angle of the source is similar at 7 mm and 1.3 cm, as well 
as to the position angle of the monopolar SiO outflow reported by Beuther et al. (2004a). 

Also the 3.6 cm continuum map shows a source that is coincident with IRAS 
18089-1732(l)a (see Figure la), however, with our resolution (~ 4") we cannot resolve it. 
From figure 1 of Beuther et al. (2004a) and our Figures 1 and 2, we associate this 3.6 cm 
continuum source with the source found previously by Beuther et al. (2002b) at the same 
wavelength. Moreover, the flux densities for both epochs are quite similar, ~1 mJy. 

Figure 3a shows the spectral energy distribution for IRAS 18089-1732(1) a from the 
centimeter to submillimeter wavelengths. It shows a "combined", two-regime spectrum 
in which the centimeter emission is dominated by a flat or slowly rising spectrum 
(q!=0.58±0.05) which can be interpreted as moderately optically thick free-free emission, 
while the millimeter and sub-millimeter emission is dominated by a component that rises 
rapidly with frequency. This component is likely to be associated with dust emission from 
a core or disk. 

Since the orientation of the position angle of IRAS 18089-1732(l)a is in the south-north 
direction (in the 1.3 cm and 7 mm observations), and its spectral energy distribution 
at centimeter wavelengths suggests that this source may be associated with a thermal 
jet or stellar wind, we interpret this source as a thermal jet that is possibly driving the 
north-south SiO molecular outflow found by Beuther et al. (2004b). Moreover, as we will 
see in our sample, these free-free emission sources of relatively flat spectrum appear to be 
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systematically associated with CH3OH maser spots (see Figure lb and discussion below). 

4.1.2. IRAS 18089-1732(l)h 

This source is reported here for the first time. IRAS 18089-1732(l)b is resolved in the 7 
mm maps. It has deconvolved dimensions of 2''7 ± 0''2 x 1''6 ± 0''2 and a P.A. = 35° ±46°. 

IRAS 18089-1732(l)b is only detected at 7 mm. Figure 3b shows its spectral energy 
distribution from the centimeter to submillimctcr wavelengths, mostly upper limits. From 
our figures and the SED this source could be interpreted as associated with dust emission. 
However, the spectrum of this source seems to flatten between 1.2 mm and 850 //m (see 
Figure 3 and Beuther et al. 2004b). Therefore, we speculate that IRAS 18089-1732(l)b 
could be also an optically thick H II region that has its turnover frequency at about 100 
GHz. 

4.2. IRAS 18182-1433 

This region contains a strong and extended (~ 13") 1.2 mm continuum source and H2O 
and Class II CH3OH masers spots (Beuther et al. 2002a, 2002c). In recent SMA millimeter 
continuum and line observations (Beuther et al. 2005, in preparation), a compact core (~ 
4") has been detected that is also associated with a molecular line forest. This millimeter 
continuum core is associated with the H2O masers. Beuther et al. (2002b) also flnd an 
extended (~ 1') C0(2-l) molecular outflow emanating from this region. 

The 3.6, 1.3 and 0.7 cm maps (see Figures 4 and 5) show that IRAS 18182-1433 
is a triple source with components separated by 2 to 10 arcscconds. Wc refer to these 
components as IRAS 18182-1433a, IRAS 18182-1433b and IRAS 18182-1433c. 
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4.2.1. IRAS 18182-1433a 

IRAS 18182-1433a is detected in our 7 mm continuum map, but not in our 3.6 cm and 
1.3 cm maps above the 4cr level of 55 and 70 ^Jy, respectively. This source is not resolved 
in the 7 mm continuum map. However, based on our synthesized beam, we set an upper 
limit of < 2 " for its angular size. 

In Figure 5 we can see that this source is associated with the water masers and the 1.2 
mm continuum core found by Beuther et al. (2002a,c). However, there is a shift of about 
0.3" between the two peaks at 7 and 1.2 mm. One possibility, is that this shift could be 
due to the fact that dust emission contribution at 1.2 mm comes from both sources, IRAS 
18182-1433a and IRAS 18182-1433b, while the source IRAS 18182-1433b was not detected 
in our 7 mm continuum map. 

Figure 6a shows the spectral energy distribution for IRAS 18182-1433a from the 
centimeter to submillimeter wavelengths. It shows a spectrum associated with a strong dust 
emission, possibly from the core and disk. 

4.2.2. IRAS 18182-1433h 

IRAS 18 182- 1433b is resolved in our 1.3 cm continuum maps. It has deconvolved 
dimensions of 1'.'7 ± 0'.'4 x 1'.'4 ± 0'.'6 and a RA. = 87° ±40°. 

IRAS 18182-1433b is not detected in our 7 mm map above a 4(T level of 0.22 mJy. The 
CH3OH Class II maser spot found by Beuther et al. (2002c) is associated with this source. 
Figure 6b shows its spectral energy distribution from the centimeter to submillimeter 
wavelengths. It shows a slowly rising (q!=0.10±0.01) spectrum. We propose that is 
produced by free-free emission from a thermal jet or stellar wind. However, this flat spectral 
index is also consistent with an optically thin H II region. 
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4.2.3. IRAS 18182-1433C 

This source is not resolved in our 1.3 cm maps. Again, based on our synthesized beam, 
we set an upper Hmit of < 2" for its angular size. IRAS 18182-1433c is not detected in our 
7 mm map above a 4(7 level of 0.22 mJy. 

Figure 6c shows its spectral energy distribution from the centimeter to submillimeter 
wavelengths. It shows a slightly negative spectrum (q;=-0.36±0.05), which is suggestive of 
non thermal radiation. Similar negative spectral indices found in sources associated with 
massive and low-mass star formation are believed to be due to synchrotron contamination 
from strong shocks (Serpens: Rodriguez et al. 1989, HH 80-81: Marti, Rodriguez, & 
Reipurth 1993, Cep A: Garay et al. 1996, W3(0H): Wilner, Reid, & Menten 1999, 
IRAS 16547-4247: Rodriguez et al. 2005 G192.16-3.82: Shepherd & Kurtz 1999). A detailed 
discussion on the presence of thermal and non thermal components in this type of objects 
is given in Garay et al. (1996). 

4.3. IRAS 18264-1152 

This region contains a strong and extended (~ 15") 1.2 mm continuum source and H2O 
and Class II CH3OH maser spots (Beuther et al. 2002a, 2002c). IRAS 18264-1152 shows 
evidence of outflowing activity. Beuther et al. (2002b) found an extended (~ 1') C0(2-l) 
molecular outflow emanating from here. 

The 7 mm and 1.3 cm maps (see Figures 7 and 8) show that this source is a triple 
source, with components separated approximately by ~2". We refer to these components as 
IRAS 18264-1152a, IRAS 18264-1152b, and IRAS 18264-1152c. We propose that they are 
possibly tracing the components of a triple stellar system. 



4.3.1. IRAS 18264-1152a 



Figure 9 shows the spectral energy distribution of this source from the centimeter to 
miUimeter wavelengths. The power law index of 2.1±0.05 suggests an optically thick H II 
region or dust emission from the core and disk. 

4.3.2. IRAS 18264-1152h 

This source is resolved in our 7 mm continuum maps. It has deconvolved dimensions 
of 3'.'8 ± 0:'5 X VIA ± OM and a P.A. = 143° ± 7°. 

Figure 9 shows its spectral energy distribution from the centimeter to millimeter 
wavelengths. The shghtly rising spectrum (spectral index of 0.27±0.06) between 3.6 and 
1.3 cm suggests a thermal jet or a partially optically thick H II region. The flux density at 
7 mm is much larger than the value extrapolated from the 3.6 and 1.3 cm observations and 
suggests that dust emission may be dominant at 7 mm. 

4.3.3. IRAS 18264-1152C 

This source is resolved in our 7 mm continuum maps. It has deconvolved dimensions 
of 2M ± 0:'5 X 0:'9 ± 0'.'7 and a P.A.= 45° ± 13°. 

Figure 9 shows its spectral energy distribution from the centimeter to millimeter 
wavelengths. At 3.6 cm, we have only an upper limit. The power law index of 2.1±0.05 
between 1.3 cm and 7 mm suggests an optically thick H II region or dust emission from a 
core and disk. 
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4.4. IRAS 18308-0841 

This region contains a strong and extended (~ 17") 1.2 mm continuum source and 
one H2O maser spot (Beuther et al. 2002a, 2002c). Beuther et al. (2002b) found strong 
^^C0(2-l) molecular emission. However the presence of bipolar outflow structure was not 
clearly estabhshed. 

The 7 mm and 1.3 cm maps (see Figures 10 and 11) show that this source is a triple 
source with components separated by 3" to 10". We refer to these components as IRAS 
18308-0841a, IRAS 18308-0841b and IRAS 18308-0841c. 

4.4.1. IRAS 18308-0841a 

This source is resolved in our 7 mm continuum maps. It has deconvolved dimensions 
of 2!'16 ± 0!'04 X 2!'0 ± 0('05 and a P.A. = 92° ± 14°. This source is also resolved in our 
1.3 cm continuum maps. It has deconvolved dimensions of 1"3 ± 0"9 x 0"7 ± 0"2 and a 
RA. = 100° ± 8°. Finally, the source is also resolved in our 3.6 cm continuum maps. It 
has deconvolved dimensions of 2"4 ± 0"1 x 2"20 ± 0"05 and a P.A. = 60° ± 16°. We find 
that while the 7 mm and 1.3 cm deconvolved dimensions are consistent, the deconvolved 
dimensions at 3.6 cm appear to be significantly larger. This may be due to the presence of 
extended emission that is resolved out or too faint to be detected at 7 mm and 1.3 cm. 

Figure 12 shows its spectral energy distribution from the centimeter to millimeter 
wavelengths. The negative spectral index of -0.42±0.07 could be suggest the presence 
of synchrotron emission. However, this spectral index could also be due to missing flux 
density in the higher angular resolution and higher frequency observations. Matching beam 
observations of the source at several frequencies are needed to discuss the spectral indices 
in a very reliable way. 
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4.4-2. IRAS 18308-0841h 

This source is resolved in our 3.6 cm continuum maps. It has deconvolved dimensions 
of 3'.'1 ± 0:'3 X 2'!b ± 0'.'2 and a P.A. = 75° ± 12°. 

Figure 12 shows its spectral energy distribution from the centimeter to millimeter 
wavelengths. Again, the negative spectral (-0.26±0.01) index could be suggest the presence 
of synchrotron emission. However, this spectral index could also be due to missing flux 
density in the higher angular resolution and higher frequency observations. Matching beam 
observations of the source at several frequencies are also needed to discuss the spectral 
indices in a very reliable way. 

4.4.3. IRAS 18308-0841C 

Figure 12 shows its spectral energy distribution from the centimeter to millimeter 
wavelengths. It was detected only at 7 mm, suggesting that dust is dominant at this 
wavelength. 

4.5. Clusters of UCHII Regions in the Vicinity of the sources IRAS 

18089-1732(1) and IRAS 18182-1433 

We report for the first time the detection of seven ultracompact HII regions in the 
vicinity of the source IRAS 18089-1732(1) and two more ultracompact HII regions in the 
source IRAS 18182-1433 at 3.6 cm (see Figure 13 and 14). Their physical parameters are 
shown in Table 5. Most of them show a comctary or unresolved spherical morphology. 
Following Van der Tak and Menten (2005), and assuming that the free- free emission is 
optically thin, that the ionized gas has a temperature of 10^ K, and that these sources have 
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distances of 3.6 and 4.5 kpc, we find that the flux of Lyman continuum photons for each 
UCHII, can be provided by an early B-type star (see Panagia 1973). 

We searched for more UCHII regions in the vicinity of the others two IRAS sources 
that were observed at 3.6 cm, however, we did not detect other cases at a Aa level of their 
respective rms noises. Therefore, in comparison with these regions, IRAS 18089-1732(1) 
and IRAS 18182-1433 can be considered rich regions of high mass star formation. 

4.6. General Discussion of the Sources 

With the present angular resolution and sensitivities of our images, it is interesting to 
note that the four sources detected at 7 mm (IRAS 18089-1732, IRAS 18182-1433, IRAS 
18264-1152, and IRAS 18308-0841) are double or triple, with no case of a single source. In 
the four regions with detected 7 mm emission, we find a total of 11 subcomponents. The 
spectral index of these subcomponents varies and suggests different interpretations for the 
nature of the sources. In Table 6 we give a tentative interpretation for the nature of the 
sources, based on our discussion. 

We note that the sources IRAS 18089-1732(l)a, IRAS 18182- 1433b, IRAS 18264-1152b 
hat exhibit slightly rising spectral indices (that are frequently associated with winds or 
thermal jets), show a correlation with methanol masers spots. Furthermore, those maser 
spots are found in the same orientation of the molecular outflows (See flgure lb, 5 and 
8). Since the methanol masers already have been suggested to be associated with shocked 
molecular gas in the outflows produced by the star (De Buizer 2003; Kurtz 2004), our data 
show this correlation, and supports the notion that the methanol masers are associated 
with ionized or shocked gas rather than disks (or dust emission). 

It is puzzling that we only detected 4 of the 10 sources observed, since they all (with the 
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exception of IRAS 18553+0414 and IRAS 18089-1732(4)) seem to have similar luminosities, 
and distances. However, these four detected continuum sources (IRAS 18089-1732, IRAS 
18182-1433, IRAS 18264-1152 and IRAS 18308-0841) are embedded in the most massive 
cores of our initial subset (see Table 1). 

Of the 11 subcomponents detected, four (IRAS 18089-1732(l)b, IRAS 18182-1433a, 
IRAS 18264-1 152b, IRAS 18308-0841c) have 7 mm emission that could be coming from the 
dust. 

5. Conclusions 

Our main conclusions are summarized below. 

1) We mapped ten regions of massive star formation at 7 mm with the VLA, detecting 
compact emission in four of them. These four sources were further observed with the VLA 
at 1.3 and 3.6 cm. 

2) Multiplicity is a common factor in the emission detected: three of the four sources 
detected are resolved into three subcomponents, while the fourth source is a double source. 

3) Four of the sources detected show 7 mm emission above the value extrapolated 
from the observations at longer wavelengths assuming jet emission. These sources are 
IRAS 18089-1732(l)b, IRAS 18182-1433a, IRAS 18264-1152b, IRAS 18308-0841c. They 
are considered good candidates of massive young stars with dust emission at 7 mm that 
could come from a disk or envelope. However, further observations at higher and lower 
frequencies are required to establish if the detected 7 mm emission is actually coming from 
dust or from optically thick ionized gas (free-free emission). 

4) Two of the sources show spectral indices consistent with ionized outflows: IRAS 
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18089-1732(l)a and IRAS 18182-1433b. IRAS 18264-1152b show, in addition to the 7 mm 
excess, a spectral index consistent with an ionized outflow between 3.6 and 1.3 cm. This 
last source may then be a combination of thermal jet and disk. 

5) Two of the sources show spectral indices consistent with optically thick HII regions: 
IRAS 18264-1 152a and c. 

6) Three of the sources show negative spectral indices that suggest the presence of 
synchrotron contamination: IRAS 18182- 1433c, and IRAS 18308-0841a and b. However, 
these spectral indices could also be due to missing flux density in the higher angular 
resolution, higher frequency observations. If conflrmed with multifrequency, matching beam 
observations, this synchrotron emission could be an indicator of strong shocks produced by 
the powerful outflows known to exist in these regions. 

8) Finally, we also find that the methanol masers in these regions are systematically 
associated with radio sources that show slightly rising spectral indices: IRAS 18089- 
1732(l)a, IRAS 18182-1433b, and IRAS 18264-1152b. This indicates that the methanol 
maser emission is associated with outflows rather than with disks. The interpretation is 
that these masers are associated with the presence of compact free-free emission. 

Our flndings are consinstent with the accretion scenario, where a dynamical collapse 
of molecular clumps, results in the formation of disks and molecular outflows, which leads 
to the formation of a massive stellar core. However, in each IRAS source we also detected 
multiple radio sources separate a few arcseconds, suggesting the formation of groups of 
massive stars. 

L.A.Z. acknowledges the Smithsonian Astrophysical Observatory for a predoctoral 
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Fig. 1. VLA continuum images of the source IRAS 18089-1732(1). A). The 3.6 cm continuum emission. The contours 

are -4, 4, 6, 8, and 10 times 40 ixJy beam~^, the rms noise of the image. B). The 1.3 cm continuum emission (contours) and 7 
mm continuum emission (grey scale). The contours are -4, 4, 6, 8, 10, 12, 14, 16 and 18 times 70 /ijy beam^^, the rms noise 
of the image. The half power contour of the synthesized beam is shown in the bottom left corner of each imagen. The 7 mm 
continuum map has a (u,v) tapering of 200 kA. The square and triangle indicate CH3OH and H2O masers positions (Walsh et 
al. 1998 and Beuther et al. 2002), respectively. The line indicates the direction of the SiO molecular outflow found by Beuther 
et al. (2004b). 
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Fig. 2 . VLA and SMA continuum images of the source IRAS 18089-1732(1). The contours shows the 7 mm continuum 

emission and the grey scale shows the SMA continuum emission at 1.2 mm (A) and 860 /im (B) (Beuther et al. 2004b, 2005). 
The half power contour of the synthesized beam is shown in the bottom left corner of each image. In the A and B images the 
contours are -4, 4, 5, 6, 7 and 8 times 0.19 mjy beam~^, the rms noise of the image. The 7 mm continuum map has a (u,v) 
tapering of 150 kA. 
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Fig. 3. Spectral energy distribution (SED) for each detected radio continuum source in IRAS 18089-1732(1) combining 

the 0.7, 1.3 and 3.6 cm VLA continuum data and the SMA continuum 1.2 and 0.7 mm data from Beuther et al. (2004a, 2005b). 
The squares are detections, the respective error bars were smaller than the squares and are not presented. The triangles with 
arrows are upper limits (4cr). The line is a least-squares power law fit (of the form Si, oc 1'°') to each spectra. In the upper 
pannel, the 1.2 mm measurement is the peak emission of this source. 
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Fig. 4. VLA continuum images of the source IRAS 18182-1433. A). The contours show the 3.6 cm continuum emission 

and the grey scale shows the emission at 7 mm. The contours are -4, 4, 5, 6, 7, 8 and 10 times 55 /iJy beam~^, the rms noise of 
the image. B) The contours show the 1.3 cm continuum emission and the grey scale shows the emission at 7 mm. The contours 
are -4, 4, 5, 6, 7, 8 and 10 times 70 fiJy beam~^, the rms noise of the image. The half power contour of the synthesized beam 
is shown in the bottom left corner of each image. The 7 mm continuum maps has a (u,v) tapering of 100 kA. The line indicates 
the North-South CO(2-l) molecular outflow found by Beuther et al. (2002b). 
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Fig. 5. VLA continuum image at 7 mm of the source IRAS 18182-1433. The half power contour of the synthesized beam 

is shown in the bottom left corner of the image. The contours are -4, 4, 5, 6, 7 and 8 times 0.22 mjy beam~^, the rms noise of 
the image. The square and triangles indicate CH3OH and H2O maser positions (Walsh et al. 1998 and Beuther et al. 2002c), 
respectively. This 7 mm continuum map has a (u,v) tapering of 100 kA. The cross indicates the 1.2 mm continuum emission 
peak found by Beuther et al. (2005), in preparation. 
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Fig. 6. Spectral energy distribution (SED) for each detected radio continuum source in IRAS 18182-1433 combinig the 

0.7, 1.3 and 3.6 cm VLA continuum data. The squares are detections, the respective error bars were smaller than the squares 
and are not presented. The triangles with arrows are upper limits (4(7). The line is a least-squares power law fit (of the form 
Sv (XI/™) to each spectra. 
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Fig. T. VLA continuum images of the source IRAS 18264-1152. A). The 3.6 cm continuum emission. The contours are 

-4, 4, 6, 8 and 10 times 40 /ijy beam~^, the rms noise of the image. B) The contours show the 1.3 cm continuum emission 
and the grey scale shows the emission at 7 mm. The contours are -4, 4, 5, 6, 7 and 8 times 60 //Jy beam~^, the rms noise 
of the image. The half power contour of the synthesized beam is shown in the bottom left corner of each image. The 7 mm 
continuum map has a (u,v) tapering of 100 kA. The line indicates the East- West CO(2-l) molecular outflow found by Beuther 
et al. (2002b). 
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Fig. 8 . VLA continuum image at 7 mm of the source IRAS 18264-1152. The half power contour of the synthesized beam 

is shown in the bottom left corner of the image. The contours are -4, 4, 5, 6, 7, 8 and 10 times 0.18 mjy beam~^, the rms 
noise of the image. The square and triangles indicate CH3OH and H2O mascrs positions (Beuther et al. 2002c), respectively. 
This 7 mm continuum map has a (u,v) tapering of 100 kA. It is interesting to note that every 7 mm radio continuum source 
has associated a H2O maser spot. 
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Fig. 9. Spectral energy distribution (SED) for each detected radio continuum source in IRAS 18264-1152 combining the 

0.7, 1.3 and 3.6 cm VLA continuum data. The squares are detections, in a few cases the respective error bars were smaller than 
the squares and are not presented. The triangles with arrows are upper limits (4cr). The line is a least-squares power law fit 
(of the form Si, oc to each spectra. 



-31 - 



-0839 00 



05 



O 
5 



U 
ui 
a 



10 



15 



20 



25 



30 



35 

18 33 34.0 



-08 39 00 



05 



U 
u 
Q 



10 



15 



20 



25 



30 



35 







33.5 





33.0 



b 



32.5 




_L 



32.0 



1833 34.0 33.5 33.0 32.5 32.0 

RIGHT ASCENSION (J2000) 

Fig. 10. VLA continuum images at 3.6 (A) and 1.3 (B) cm of the source IRAS 18308-0841. The half power contour of 

the synthesized beam is shown in the bottom left corner of each image. A). The contours are -4, 4, 6, 8, 10, 12, 14, 16, 18, 
20 and 30 times 0.16 mjy beam~^, the rms noise of the image. B). The contours are -4, 4, 6, 8, 10, 12, 14, 16, 18, 20 and 30 
times 0.17 mJy beam~^, the rms noise of the image. 



-32- 



-0839 00 



05 



g 

5 



u 

lU 

□ 



10 



15 - 



20 



A 



25 



30 



35 



1833 34.0 



_L 



_L 



33.5 33.0 32.5 

RIGHT ASCENSION (J2000) 



32.0 



Fig. 11. VLA continuum image at 7 mm of the source IRAS 18308-0841. The half power contour of the synthesized 

beam is shown in the bottom left corner of the image. The contours are -4, 4, 5, 6, 7, 8 and 10 times 0.20 mjy beam~^, the 
rms noise of the image. The triangle indicates the H2O maser position (Beuther et al. 2002c). This continuum map has a (u,v) 
tapering of 200 kA. The dot indicates the MAMBO 1.2 mm continuum emission peaJc found by Beuther et al. (2002a) 
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Fig. 12. Spectral energy distribution (SED) for each detected radio continuum source in IRAS 18308-0841 combining 

the 0.7, 1.3 and 3.6 cm VLA continuum data. The squares are detections, in most of the cases the respective error bars were 
smaller than the squares and are not presented. The triangles with arrows are upper limits (4(t). The line is a least-squares 
power law fit (of the form Sn oc i^") to each spectra. 
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Fig. 13. VLA continuum image at 3.6 cm towards the source IRAS 18089-1732(1). This map was made with ROBUST=5 

(equivalent to natural weighting). The half power contour of the synthesized beam is shown in the bottom left corner of the 
image. The contours are -3, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 30 and 60 times 40 /iJy beam~^, the rms noise of the image. The 
sources G12.879+0.496, G12.891+0.495, G12.886+0.487, G12.918+0.487, G12.918+0.487, G12.910+0.477, and G12.944+0.488 
are first reported here. This map was not corrected by primary beam response. The box indicates the approximate size of the 
figure la. 



-35- 



-14 28 30 
29 00 
30 



o 30 00 
o 

CM 



< 



O 
UJ 
Q 



30 



31 00 



30 - 



32 00 



30 



33 00 a 



G1 6.61 4-0.01 8 



IRAS 18182-1433 



1821 10 



05 



00 



G1 6.560+0.001 




20 55 



RIGHT ASCENSION (J2000) 



Fig. 14 . VLA continuum image at 3.6 cm towards the source IRAS 181821433. This map was made with ROBUST=5 

(equivalent to natural weighting) . The half power contour of the synthesized beam is shown in the bottom left corner of the 
image. The contours are -3, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 30 and 60 times 3.5 /iJy beam~^, the rms noise of the image. The 
sources G16. 560+0. 001 and G16. 614-0. 018 are first reported here. This map was not corrected by primary beam response. The 
box indicates the approximate size of the figure 4a. 
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Table 1. Physical Parameters of the 10 Selected 1.2 mm Continuum Sources' 



Physical Parameters 



Coordinates Distance^ Luminosity Core Mass'^ 

Source n-2non '^2non [Kpc] [L-] [M : ] 



IRAS 


18089-1732(1)^ 


18 


11 51.3 


-17 31 29 


3.6 


3.2 


X 


10^ 


1200 


IRAS 


18089-1732(4) 


18 


12 30.4 


-17 32 59 


3.6 








200 


IRAS 


18090-1832 


18 


12 01.9 


-18 31 56 


6.6 


1.3 


X 


10^ 


1000 


IRAS 


18182-1433 


18 


21 07.9 


-14 31 53 


4.5 


2.0 


X 


10^ 


1500 


IRAS 


18264-1152 


18 


29 14.3 


-11 50 26 


3.5 


1.0 


X 


10^ 


2200 


IRAS 


18290-0924 


18 


31 44.8 


-09 22 09 


5.3 


2.5 


X 


10^ 


800 


IRAS 


18308-0841 


18 


33 31.9 


-08 39 17 


4.9 


1.6 


X 


10^ 


1300 


IRAS 


18521+0134 


18 


54 40.8 


01 38 02 


5.0 


1.3 


X 


10^ 


500 


IRAS 


18553+0414 


18 


57 52.9 


04 18 06 


0.6 


2.5 


X 


102 


15 


IRAS 


19012+0536 


19 


03 45.1 


05 40 40 


4.6 


1.6 


X 


10^ 


500 



Note. — Units of right ascension are hours, minutes, and seconds, and units of decHnation are 
degrees, arcminutes, and arcseconds. 

^These parameters were obtained from Sridharan et al. (2002). 

'^This distance Usted is the nearest distance, and the luminosity and core mass were calculated with 
this distance. 

'^The masses listed here are a factor of two less than calculated by Beuther et al. (2002a). This 
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difference is due to Beuther et al. (2002)a use a sightly different procedure to calculate the total mass 
of the cores that the developed by Hildebrand 1983. See the erratum for this paper. 

*^This source is previously known as IRAS 18089-1732, however we have added a numbering (1) at 
the end of the name to distinguish between this component and component IRAS 18089-1732(4) 
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Table 2. Parameters of the VLA 0.7 cm Continuum Observations. 



Bootstrapped Rms Synthesized Beam 







Phase Center 


Phase 


Flux Density 


Noise 


Size 


P.A. 


Source 




0^2000 


^2000 


Calibrator 


[Jy] 


[mJy] 


[arcsec] 


[deg.] 


18089-1732(1) 


18 


11 51.3 


-17 31 29 


1820-254 


0.479+0.008 


0.18 


1.63x1.08 


-29 


18089-1732(4) 


18 


12 30.4 


-17 32 59 


1820-254 


0.479+0.008 


0.16 


2.11x1.47 


-33 


18090-1832 


18 


12 01.9 


-18 31 56 


1820-254 


0.479+0.008 


0.18 


1.99x1.49 


-30 


18182-1433 


18 


21 07.9 


-14 31 53 


1832-105 


0.656+0.006 


0.19 


1.90x1.59 


-30 


18264-1152 


18 


29 14.3 


-11 50 26 


1832-105 


0.656+0.006 


0.18 


1.84x1.62 


-28 


18290-0924 


18 


31 44.8 


-09 22 09 


1832-105 


0.656+0.006 


0.16 


1.37x1.14 


-10 


18308-0841 


18 


33 31.9 


-08 39 17 


1832-105 


0.656+0.006 


0.22 


1.77x1.64 


-14 


18521+0134 


18 


54 40.8 


01 38 02 


1851+005 


0.79+0.01 


0.13 


1.35x1.18 


-9 


18553+0414 


18 


57 52.9 


04 18 06 


1851+005 


0.79+0.01 


0.13 


1.38x1.20 


-22 


19012+0536 


19 


03 45.1 


05 40 40 


1851+005 


0.79+0.01 


0.12 


1.36x1.19 


-30 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcmin- 
utes, and arcseconds. 
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Table 3. Parameters of the 1.3 cm VLA Continuum Observations. 











Bootstrapped 


Rms 


Synthesized Beam 




Phase Center 


Phase 


Flux Density 


Noise 


Size 


P.A. 


Source 


Ct2000 


(^2000 


Cahbrator 


[Jy] 




[arcsec] 


[deg.] 


18089-1732(1) 


18 11 51.3 


-17 31 29 


1820-254 


0.626±0.003 


70 


1.73x0.94 


15 


18182-1433 


18 21 07.9 


-14 31 53 


1832-105 


0.981±0.006 


70 


2.35x1.63 


28 


18264-1152 


18 29 14.3 


-11 50 26 


1832-105 


0.981±0.006 


60 


2.38x1.63 


31 


18308-0841 


18 33 31.9 


-08 39 17 


1832-105 


0.981±0.006 


170 


2.44x1.63 


36 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcmin- 
utes, and arcseconds. 
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Table 4. Parameters of the VLA 3.6 cm Continuum Observations. 











Bootstrapped 


Rms 


Synthesized Beam 




Phase Center 


Phase 


Flux Density 


Noise 


Size 


P.A. 


Source 


Ct2000 


(^2000 


Cahbrator 


[Jy] 




[arcsec] 


[deg.] 


18089-1732(1) 


18 11 51.3 


-17 31 29 


1820-254 


0.771±0.002 


40 


4.70x2.75 


17 


18182-1433 


18 21 07.9 


-14 31 53 


1832-105 


1.415±0.004 


55 


4.75x2.62 


24 


18264-1152 


18 29 14.3 


-11 50 26 


1832-105 


1.415±0.004 


40 


4.90x2.66 


31 


18308-0841 


18 33 31.9 


-08 39 17 


1832-105 


1.415±0.004 


160 


4.40x2.68 


28 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcmin- 
utes, and arcseconds. 
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Table 5. Physical Parameters of the HII Regions Detected with the VLA at 3.6 cm". 



Physical Parameters 



Source 


02000 


^2000 


Deconvolved Size 


P.A. 


Flux Density 
[mJy] 


logCNi)*- 
s-i 


G12.879+0.496 


18 11 48.834 


-17 31 49.07 


T.'4 ± 0'.'5 


X 5'.'2 ± o;'5 


75° ± 10° 


6.0+0.1 


45.7 


G12.891+0.495 


18 11 50.422 


-17 31 14.12 


3'.'0 ± 0^3 


X 2;'3 ± 0{'6 


115° ±68° 


3.0+0.1 


45.5 


G12.886+0.487 


18 11 51.510 


-17 31 40.57 


6:'4 ± 1^4 


X 3;'7 ± 1{'5 


58° ± 25° 


1.0+0.1 


45.0 


G12.883+0.480 


18 11 52.680 


-17 32 04.92 


9'.'0 ± O'.'S 


X 8;'o ± o;'7 


16° ± 50° 


6.0+0.1 


45.7 


G12. 918+0.487 


18 11 55.360 


-17 30 01.92 


e'.'o ± of'5 


X 5;'3 ± 0"4 


44° ± 26° 


6.0+0.1 


45.7 


G12. 910+0.477 


18 11 56.751 


-17 30 44.05 


9{'6 ± I'.' I 


X 7{'0 ± 0'.'9 


42° ± 20° 


5.0+0.1 


45.6 


G12.944+0.488 


18 11 58.314 


-17 28 38.78 


3'.'6 ± o;'6 


X i;'6 ± o;'7 


43° ± 13° 


1.0+0.1 


45.1 


G16.560+0.001 


18 20 54.977 


-14 31 39.84 








2.5+0.1 


45.5 


G16.614^0.018 


18 21 05.386 


-14 29 20.45 








0.9+0.1 


45.1 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, 
arcminutes, and arcseconds. 

''These parameters were obtained from a map with the primary beam response correction applied. 

''The flux of Lyman continuum photons was obtained assuming distances of 3.6 kpc and 4.5 kpc. 



Table 6. Flux Densities at 0.7, 1.3 and 3.6 cm of the VLA Continuum Sources. 



Flux Density 
3.6 cm 1.3 cm 0.7 cm 



Source 02000 52000 [™Jy] [mJy] [mJy] Nature 



18089- 


-1732(l)a 


18 


11 


51.451 


-17 


31 


28.85 


1.10±0.05 


1.60±0.04 


3.00±0.30 


Thermal Jet -|- Core 


18089- 


-1732(l)b 


18 


11 


51.571 


-17 


31 


27.89 


< 0.20 


< 0.,30 


2.10±0.,30 


H II Region? 


18182- 


-1433a 


18 


21 


09.014 


-14 


31 


47.62 


< 0.20 


< 0.35 


2.26±0.30 


Core+Disk 


18182- 


-1433b 


18 


21 


09.129 


-14 


31 


48.59 


0.60±0.05 


0.80±0.05 


0.72±0.10 


Thermal Jet 


18182- 


-1433c 


18 


21 


09.469 


-14 


31 


56.93 


1.65±0.03 


1.10±0.05 


<0.80 


Strong Shocks? 


18264- 


-1152a 


18 


29 


14.203 


-11 


50 


23.80 


< 0.20 


0.34±0.10 


1.30±0.20 


H II Region? 


18264- 


-1152b 


18 


29 


14.356 


-11 


50 


22.50 


0.80±0.02 


1.06±0.07 


2.70±0.20 


Thermal Jet + Core 


18264- 


-1152c 


18 


29 


14.406 


-11 


50 


24.77 


< 0.20 


0.44±0.10 


1.70±0.30 


H II Region? 


18308- 


-0841a 


18 


33 


32.891 


-08 


39 


18.47 


11.40±0.05 


10.30±0.20 


5.38±0.50 


Strong Shocks? 


18308- 


-0841b 


18 


33 


32.904 


-08 


39 


07.33 


3.21±0.05 


2.46±0.10 


<0.80 


Strong Shocks? 


18308- 


■0841c 


18 


33 


33.120 


-08 


39 


07.53 


< 0.60 


< 0.75 


2.20±0.40 


Core? 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are 
degrees, arcminutes, and arcseconds. 



